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SUMMARY 

The supersonic flutter of wide panels on discrete flexible supports is investigated 
for three different panel- support configurations. The first study examines the effect of 
support stiffness on the flutter behavior of a specific five-support configuration with 
leading- and trailing-edge overhangs; this configuration was recently under consideration 
for a heat shield on an entry body. The second study investigates the effect of support 
stiffness on flutter of panels with various numbers of equally spaced supports. The third 
study examines the effect of center- support location on the flutter of panels with three 
supports. Results are presented in nondimens ional form. The analysis is based on wide- 
plate structural theory and Ackeret aerodynamics. Finite differences are employed to 
obtain solutions for flutter pressure. A computer program based on this analysis and 
including a direct solution technique is presented. The program can be used to find the 
flutter pressure of wide panels of variable thickness supported by any number of flexible 
supports. 


INTRODUCTION 


Panels held in place with rows of standoff supports, as shown in figure 1, have been 
considered for thermal protection systems of vehicles designed to withstand reentry from 
space. Static design of such systems generally leads to panels of minimum gage that are 
held in place by rather flexible supporting structure. However, dynamic considerations, 
particularly panel flutter, must also be satisfied and may dictate panel gage as well as 
dimensions of the supporting structure. Panel flutter may also be a design consideration 
in aircraft and other vehicles in which the skin panels often have lines of discrete supports. 

A general discussion of flutter of panels is presented in reference 1, Results for 
the flutter of wide panels on nondeflecting supports are presented in reference 2; however, 
no analytical results are available for the flutter of panels on multiple discrete flexible 
supports. Wide-plate structural theory coupled with Ackeret aerodynamic theory can be 
used to obtain meaningful estimates of supersonic flutter speeds (dynamic pressure) for 
many configurations of interest with a minimum amount of computation. Additionally, 



discrete flexible supports can be readily accounted for through the method of writing 
equilibrium in terms of finite differences. Thus, in order to provide flutter predictions 
for wide panels on multiple flexible supports, equilibrium equations were developed in 
conjunction with this combination of theories (see "Method of Analysis" and appendix A), 
and a computer program was written which employs a direct solution technique rather 
than the usual modal approach. A listing of this program, which can be used for a wide 
class of flutter problems, is presented in appendix B. 

This paper presents three parameter studies obtained by use of the computer pro- 
gram. The first study examines the effect of support stiffness on flutter behavior of pan- 
els in the configuration of figure 1. In the second study, the effect of support stiffness on 
flutter of panels with various numbers of equally spaced supports is examined. In the 
third study, the effect of the location of the center support of a three- support panel is 
examined. 

SYMBOLS, 


a length of panel 

C support rotational stiffness 

c speed of sound 

D plate bending stiffness 

j,j n ,n integers 

k support stiffness per unit width 

M Mach number 

M X ,M X bending moments per unit width (see eqs. (A5) and (A2)) 
m mass per unit area 

N number of rows of supports 

q dynamic pressure 
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t 


time 


w,W deflections of panel 

x coordinate in lengthwise direction 

(3 = \Jm 2 - 1 


5jj n Kronecker delta 

6(x-x n ) Dirac delta function 

e spacing of nodal points 


p air density 

co frequency 

Primes denote derivatives with respect to x. 


METHOD OF ANALYSIS 


The structure considered is a wide plate with the length interrupted by lines of 
spring supports (fig. 1). The plate is in a supersonic airstream with airflow perpen- 
dicular to the lines of support. The loading considered includes inertia and aerodynamic 
loading appropriate for flutter at Mach numbers of about 1.4 and higher. (See ref. 2.) In 
this section of the paper, terms in the differential equation of equilibrium will be identi- 
fied, and the method of solution will be indicated. 

The differential equation of equilibrium of forces for a wide panel on N supports is 


3^Mx 


N 

2 

n=l 


( ~ 9 2 w\ R/ x 9^w 2q 9w 9w _ 

kw - C — r- 6(x-x n ) + m — - + — — + pc — = 0 


8X 4 


'at 2 0 3x 


9t 


( 1 ) 


where M x = D — The first term of equation (1) represents the restoring forces due 
9x z 

to the bending stiffness D of the panel. The terms in the summation represent the 
restoring forces due to the N rows of springs, each having an extensional stiffness k 
and a rotational stiffness C. The next term represents the inertia force. The other 
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terms are the aerodynamic lift associated with the local angle of attack — and the 
time rate of change in deflection — . Both these terms come from quasi- steady Ackeret 

at 

theory (or piston theory with compressibility correction). Aerodynamic damping effects, 
which are represented by the last term, are expected to be small for most cases (ref. 3). 

For flutter, periodic motion may be assumed and equation (1) is thus converted to 
an ordinary differential equation. This equation together with the boundary conditions 
is put in finite-difference form by defining nodal points and by using conventional central 
differences. The stiffnesses and the mass were considered to be functions of x, which 
introduces no important complications with the use of finite differences. Thus, solution, 
of variable-cross-section beams supported by various springs is permitted. The details 
of the derivation are presented in appendix A, which also includes the method of solution 
and in appendix B, the computer program that was used in obtaining the results of the 
parameter studies in this paper. The method of solution used is direct, not modal, and 
requires the evaluation of a second-order determinant independent of the number of 
unknowns in the problem; the method then makes use of determinant plotting to determine 
the flutter pressure. 


PROBLEMS STUDIED 

All the configurations studied are wide panels and the airflow is along the length of 
the panel. 


Special Five -Support Configuration 

The first problem considered is the determination of the flutter pressure of a panel 
with five equally spaced supports and with overhangs like the panel of figure 1. The lead- 
ing overhang is 10.7 percent of the support spacing and the trailing overhang is two-thirds 
of the leading overhang. The supports are flexible in extension but free to rotate. Flutter 
pressure is determined for a full range of support stiffness. 

Panels On Equally Spaced Supports 

The second problem considered is the determination of the flutter pressure of pan- 
els with two, three, four, or five equally spaced supports, where the outer supports are at 
the leading and trailing edges of the panel. Again, the supports are flexible in extension 
but free to rotate and flutter pressure is determined for a full range of support stiffness. 
Comparison of the results for the five- support case with the results for the special con- 
figuration indicates the effect of the overhangs on the flutter pressure. 
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Panels On Unequally Spaced Supports 

A third set of problems considered is the flutter of panels with simple supports at 
the ends and with an additional simple support at various locations. Here the flutter 
pressure is determined as a function of the location of the additional support. 

RESULTS AND DISCUSSION 

All the problems were analyzed in terms of nondimensional parameters so that 
they apply to a wide class of dimensions within the specified conf iguration. 

Special Five-Support Configuration 

Vibration frequencies were determined for the configuration of figure 1 and the 
first 11 natural frequencies are plotted in figure 2 as a function of support stiffness. 

The natural frequencies group together in fours, and pairs within each group appear 
to touch. As the support stiffness increases, the frequencies increase until they reach 
values appropriate for a panel over simple supports. Vibration frequencies were deter- 
mined by setting the air dynamic pressure and damping equal to zero and by searching 
for values of the frequency that satisfy the equations of equilibrium. Thirteen stations 
between supports were deemed sufficient and were used in the calculations for vibration 
frequency and flutter pressure. Flutter pressures were determined, as indicated in 
appendix A, by examining the behavior of the dynamic pressure with frequency. Typical 
dynamic-pressure — frequency curves are shown in figure 3 for selected values of support 
stiffness; considerable overlap of curves is indicated. This behavior was noted and dis- 
cussed in reference 2 for a panel on inflexible multiple supports. 

The dynamic pressure at flutter is given in figure 4 as a function of support stiff- 
ness. The flutter pressure is zero for zero support stiffness and generally increases 
with support stiffness, depending upon the aerodynamic damping (altitude). The curve 
for zero damping in figure 4 is based on steady-state aerodynamics, and the flutter 
dynamic pressure according to this curve drops sharply in the neighborhood of the stiff- 
ness at which the first and second natural frequencies (fig. 2) appear to be equal. In the 
neighborhood of this sharp drop, aerodynamic damping becomes important and modest 
amounts of damping wash out the dropoff. Presumably other stiffnesses for which the 
natural frequencies appear to be equal may have associated with them other sharp drops 
in flutter dynamic pressure for zero damping. It is to be expected that such drops would 
again be washed out by damping effects, and accordingly, they are not of concern. 

Flutter modes were determined for the same four support- stiffness values consid- 
ered in figure 3 and are presented in figure 5, where they exhibit remarkable change in 
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character with change in support stiffness. Previous experience has indicated that 
flutter deflections are prominent in the trailing-edge bays. However, it can be seen 
from figure 5 that when the supports are weak, deflections are more prominent in the 
leading-edge bays. 


Panels on Equally Spaced Supports 

For steady- state theory the dynamic pressure at flutter was determined for various 
values of support stiffness and plotted in figure 6 for panels on two, three, four, or five 
equally spaced supports. For this case the outer supports are at the ends of the panel. 
Again 13 stations were placed between supports. The curves generally have the same 
shape as that discussed earlier for the special five- support configuration. The panel 
with two supports, one at each end of the panel, is an exception in that it exhibits no 
sharp dropoff. At higher values of support stiffness, the flutter dynamic pressure for 
each case stops increasing and approaches the corresponding value given in reference 2 
for panels over nondeflecting (simple) supports. To compare the results for the panel 
on five equally spaced supports of figure 6 with the results for the special configuration 
of figure 4, the length of the panel of figure 4 should be changed by the factor 112/117 
to account for the overhangs. The most interesting difference noted is that the dropoff 
extends over a larger range of support stiffness when there is no overhang. Examina- 
tion of the results shows that the dropoff occurs when the support- stiffness parameter 
ka^/D(N-l)3 is about 100; this dropoff always occurred in the neighborhood of two equal 
eigenvalues. The maximum flutter pressure attainable occurs when the value of the 
support parameter ka^/D(N-l)3 is greater than 2000. 

Panels on Unequally Spaced Supports 

The dynamic pressure at flutter was determined for panels with simply supported 
ends and with an additional simple support at various locations. Flutter dynamic pres- 
sure is plotted in figure 7 as a function of the location of the support. The highest flutter 
pressure occurs when the support is at the center of the panel, but it drops off sharply, 
picks up again, then drops continuously as the support is moved toward either end of the 
panel. It should be mentioned that the flutter pressures presented in figure 7 were 
obtained from a finite-difference representation with 31 nodal points, and no attempt was 
made in this problem to find special behavior that might appear as the interval between 
nodal points is decreased. 


CONCLUDING REMARKS 

A method of analysis has been developed for the flutter behavior at supersonic 
speeds of wide panels with discrete flexible supports. Equations developed from wide- 



plate structural theory coupled with Ackeret aerodynamic theory are solved directly by 
means of finite differences. The listing is given of a computer program which is based 
on this approach and which is applicable to flutter of panels with thickness and mass 
variations along the length of the panel; panel support is represented by discrete exten- 
sional and rotational springs across the panel at various locations along the length of the 
panel. 

Results are presented for three parameter studies. The first and second studies 
are concerned with the effect of support stiffness on panels with equally spaced supports. 
The panel flutter pressure generally increases with increasing support stiffness up to a 
limiting value, which occurs at values of the support- stiffness parameter ka^/D(N-l)3 
above 2000, where N is the number of supports, k is the support stiffness per unit 
width, a is the panel length, and D is the plate bending stiffness. At values of support 
stiffness corresponding to two equal eigenvalues, ka3/D(N-l)3 of about 100, however, 
the flutter pressure decreases to a value that is highly dependent on the damping. For 
low damping (high altitude) this value may be so low that panel flutter poses a severe 
design constraint; thus, for low damping, such values of support stiffness should be 
avoided. The third study is concerned with an end-supported panel with an additional 
off-center support; it was found that the calculated flutter pressure varies in a nonuni- 
form way with change in the location of the off-center support. 

The method of analysis presented herein was found to be adequate for the problems 
treated. This demonstrates that direct techniques may be used for flutter analysis - 
that modal techniques are not necessarily required. The use of the computer program 
for these flutter problems requires considerable interface time on the part of the analyst. 
It is difficult to automate a computer program for flutter of a panel over multiple supports 
because of the multiplicity of eigenvalues in the critical range and because of the narrow 
ranges of support stiffness where damping is important. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., December 6, 1973. 
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APPENDIX A 


ANALYSIS 

Difference equations will be derived to replace the differential equation of equilib- 
rium and the boundary conditions, a method of solution will be discussed, and in appen- 
dix B a computer program will be presented which uses the method of solution to solve 
the difference equations. 


Difference Equations 

For simplicity the aerodynamic damping terms and rotational spring term are 
omitted from equation (1), which then becomes 


3 2 M x y' 

— + 2, kw6(x-x n ) 

8x n=l 


+ m 


3 2 w 2q 9w _ n 
8t 2 P 9 x 


where 


(Al) 


m x =d2% 

dx l 


(A2) 


Periodic motion is assumed as 


w = W(x) sin cot 


(A3) 


so that equation (Al) becomes the following differential equation: 

N " 

My + Y kW6(x-x n ) - mco% + — W’ = 0 
Zj 8 

n=l 


(A4) 


with 


M x = DW" 


(A5) 


Boundary conditions considered at x = 0 or a are 


Free: 


M x = M x = 0 


(A6) 
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Simply supported: 

W = M x = 0 (A7) 

Clamped: 

W = W' = 0 (A 8 ) 

In addition, a spring could be located at x = 0 or a. 

Using nodal points with spacing e and conventional central differences gives 
equations (A4) and (A5) in finite-difference form as follows: 

N 

( M x,j +1 - 2 M x,j + M x,j-l)p + £ k n W j J 6 jj n ' m j w2w j 

n -1 

* e j |( w i+i - w 3-i)r = 0 (A9> 

with 

M x,j' D i( w i + i- 2W i +w i-i)i < A10 > 

where j n is the value of j corresponding to the location of the nth spring and Ej 
is a coefficient (usually equal to 1 ) inserted for convenience in satisfying boundary condi- 
tions yet retaining the form of equilibrium equations (A9) and (A 10) as required by the 
computer program. Besides retaining the form of equations (A9) and (A10), the boundary 
conditions are set up so that W .3 is the lowest subscript, nonzero value of W; and the 
first equation considered, therefore, is the one specifying equilibrium about the point 
corresponding to j = 3. Note that D, k, and m are now identified with subscripts, 
so that stiffness (or thickness) and mass variation may be permitted and the various 
spring supports may have different stiffnesses as well as different locations. 

For a spring- supported edge at x = 0 or j = 4, the boundary conditions from 
equation (A 6 ), including a term for the spring, become 


0 

II 

sf 

(All) 


(A 12) 
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After substitution from equation (A12), equation (A9) representing equilibrium at 
j = 4 becomes 

(2M Xj5 - 2M x>4 )1 + 2k 0 W 4 i - M 4 o;2w 4 + E 4 |(w 5 - W 3 )I * 0 (A13) 

The form of equations (A9) and (A10) may be retained and yet the boundary conditions of 
equations (All) and (A 13) are satisfied if Wg, D 4 , m 4 , and E 4 are the first nonzero 
values of these parameters according to subscript, D 4 and m 4 are one -half their 
corresponding interior value, and E 4 = -i. The initial equation to be considered corre- 
sponds to j = 3. In the interior, the values of Ej are always unity. For a free edge 
at j = 4, simply set kg = 0. 

For a free edge at x = 0, halfway between j = 3 and ] = 4 die boundary condi- 
tions from equation (A 6 ) become 


m x,3 + m x,4 

~ - 


0 


(A14) 


M. 


- M 


x >3 x,4 _ 


(A15) 


or M x 3 = M x 4 = 0 . For this condition, the form of equations (A9) and (A10) is retained 
by starting with the same values as before except that D 4 , m 4 , and E 4 are the same 
as their corresponding interior values. 

For a simply supported edge at x = 0 or j = 2, the boundary conditions from 
equation (A7) are 

w 2 = Mx >2 = 0 (A16) 

The form of equations (A9) and (A10) is retained for these conditions if the first nonzero 
values of the parameters according to subscript have subscript 3 and the initial equation 
corresponds to j = 3. 

For a clamped edge at x = 0 or j = 2, the boundary conditions from equation (A 8 ) 
are 


W 2 = 0 (A17) 

W 3 - W x = 0 (A18) 
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Thus, 

2D 2 

M x ,2 = -o- W 3 (A19) 

To retain the form of equations (A9) and (A10), the first nonzero values of the parameters 
according to subscript are W 3 , D 2 , m 3 , and E 3 , and the initial equation corresponds 

to j = 3. The value taken for D 2 must be twice its corresponding interior value. 

Similar treatment must be given to boundary conditions at x = a. 

Method of Solution 

The difference equations just presented are linear and homogeneous (right-hand 
side = 0 ) and each has a maximum of five unknown deflections Wj ; the values of j 
that appear in each equation are adjacent, that is, j + 2, j + 1, j, j - 1, and j - 2. The 
equation can be arranged so that the first and last equations have only three unknowns and 
the second and next to last equations have only four unknowns. When written in matrix 
form, they can be identified as being banded, with a bandwidth of 5. A modal solution would 
start out by setting q = 0 and then solving for the modes and frequencies for which the 
determinant of the coefficients vanishes. The modes and frequencies would then be used 
in a separate flutter analysis. A direct solution would usually be obtained by finding values 
of q and co for which the determinant of the coefficients vanishes. The present method 
of solution is direct, but instead of working with the determinant of the coefficients of the 
banded matrix, a marching procedure is used. 

If the first two unknowns (W3 and W4) are assumed, the first equation will deter- 
mine the third unknown, the second equation will determine the fourth unknown, and so 
forth, and all the unknowns would be determined without using the last two equations. The 
present method of solution for given values of q and co assumes two linearly indepen- 
dent sets of values for the first and second unknowns, and thus, two sets of preliminary 
solutions are obtained on the basis of these assumptions. Then a linear combination of 
the preliminary solutions is substituted into the last two equations and the determinant 
of the coefficients of these two equations is found. A solution curve is obtained by finding 
combinations of q and u> for which this determinant is zero. 

The standard coalescence approach for undamped flutter was used in this method 
of solution in that extremum values of q along the q-co solution curves correspond to 
flutter values (fig. 3), with the lowest of such values determining the flutter pressure of 
interest. For the damped case, the solution points where the real and imaginary parts 
of the determinant vanish simultaneously correspond to flutter values. 
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COMPUTER PROGRAM 

The computer program based on the equations just derived in appendix A applies 
to wide panels of any prescribed stiffness distribution, any number and location of dis- 
crete flexible supports, any number and location of discrete rotational supports, and 
any mass distribution. 


Because of the intersecting solution curves that occur for the panels of interest, 
it is difficult to automate the solution for the flutter pressure. However, the computer 
program has options to determine q for a set of co, co for a set of q, the value of 
the determinant for each o?-q combination, and for this option, an estimate of o) 
when there is a determinant crossing at a given value of q. 

P is j G K AM SuPRFL X ( INPUT , OUTPUT , T AP£5 = INPUT ? TAPE 6-OUT PUT ) 
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* * * 
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C SUPPORT ROTATIONAL ST IFPNcSS rC* A/ C*EL 
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NAMFlI ST/PLUT gR/CMEGAi , DELCM , QMF I N , Q , CELQ , QF I N 
NAMclIST/STIFF/E I, £K,C,LM, Ew 
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A tiSCI SA { X3 tX2,Xl«Y2,Y2«Yi I--X3+ ( ( Y1*Y3-Y 2*Y3 I *( X3—X1 )*( X2-X3) >/(Yl* 
1Y2*<X1-X2) + U*Y3*<.X3-Xil*Y2*Y3*.{X2-X3) ) 

I RtiAD *5, 50 Ul ) CASE 
IF<£JF,5.I 4,5 

4 PRINT 6997 

6997 FORMAT*//* S7.UPP.fcO ON EOF, 5 *) 

STOP 

5 WRITE ( o » 5 000 1 CASE 
Re AD OPTION 

READ FLUTCR 
WRITE < 6 , OPT I ON) 

WKlTc (6, f OUTER ) 

N:u£c=N.DS +4 
READ STIFF 
U ME GA 4 . - 0 Mfc G A i / EL * * 2 
OMEGA = DM 6GA1 
DLLGM=Dfc LQM/EL **2 
GMF I N= UMF I N/ LL **2 
w=Q/EL**5/2. 

DEL Q = Ot LQ/ EL** 3/ 2 * 
wFIN-wFIiM/EL**^/2. 

5 0 DU FORM ATC////15H INPUT FOR CASE/// 1X8A1Q ) 

5001 FORMAT* 6A10> 

WRI Te I6,5o09)i I ,E I ( 1 1 ,EMU i,EKU ) ,C( I) ,EQ(I ) , I~i, NDEE ) 

5009 F OR M ATI/ 5 n 0 I 6X5 HE i (I ULX5H6M( I U*X5FbK( I) 12X4HC (I) 11 X5HEQ l I) 
i / i I 5 » 5 2 1 6 • 8 ) ) 

Ku UN T = 0 
W K I T l; ( 6 , 5 u 1 0 ) 

5 0 4 . 0 FORMAT* /9x5HUMEGA16XxHi a ii9X^hDcT^4Xl4HN0* ITERATIONS/) 

IFliJPT. 20 • I ) GO TO luu 
1F( I JPT. ;Q.3) GO TO 3uO 

C DcTLRHlNL g FOR lACH GMlU A 

30 KUUNT- u 
HuL 0 = 0 
DO 4 J 1=1,2 
I F ( 1 • Gl • 2 ) U = 0 F D.L L W 

cal l DL T cRM ( Dm E G A , g , t M-, c I f NU S f D £ T , C , E K , W , E Q ) 

X ( l ) -O.fcT 

Y l I ) = o 

40 CJmTINOc: 

50 Cl SSA = AbSClSA( Y ( 3 ) ,YL2), *U> ,a(j) , X(2) ,XU) ) 
g=CISSA 

I F { K U UN T • v W • o ) XC = X i 1 ) 

KOUNT = kOUnT + 1 

CALL UL T l RH ( OMEGA, W, tM , L I , ADS, GET ,C , E K , W , EQ ) 

A(x ) =DLT 

Y i i = 0 

1 F { A 6 S ( X ( 1 ) .) / A 6 S ( a 0 i « L c # a * 6 ~ 5 ) bO TO 3 0 
I F t NO UN 1 • GH « * 0 ) Gu TO 7 u 
U c L A b S ( ( Q— HO L D ) / 1 0 • I 
DO 60 I = L, 3 
g=Q+OLLC 

CALL UcTl;RFVlGM'cGA,Q, EM, Elf NOS, DE T , C , EK , W , EQ ) 

a ( I )-DLT 
YU ) -Q 

60 CONTI NUc 
GU Tu 50 
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70 0=0M£gA*EL**2 

QQ;=Q*£l** 5 * 2 . 

wRITu (6,2501 C, CGfDfcT.KCLNT 
XRI Tfc ( G , 2oU ) 

G u TO j. 

80 P<=0MfcGA*EL **2 
0u = Q iis EL**j; > ! e 2» 

i«RlT£(b,25G)0,GQ»0fcT,KCUNT 
UMEGA=0 MEgA+ 0ELCM 
IF t l)M EGA . GT . UMF I N ) GO TC 1 
GU TO v*u 

C DbTeRMINfc CLT FUR cACH (GmEgA,Q) 


100 L =0 
1C1 L = L + _ 

CALL DeT SRM ( CMEG A , Q , EM » e I , i\Og , DET ,C » EK , W , EQ ) 

OMcG( L)=GMEGA*£L**2 
Oe TT ( L ) =DfcT 
QQ=U*F.:l**:>*2. 

WRITE! o, 250 ) UMfcG ( L ) »QU»DLTT (L) »KOUNT 
IF ( L. LT. 3) GU TO lu2 

C IDENTIFY CHANGE IN SIGN 

Ql = OETT(U*DbTT(L-l) 

D2=DETT (L) *OLTT ( L-2) 

1 F ( 01 « GE ® o « A NO • C c • G E « 0 ) GC TL a 0 2 

C INTERPuLATt 

i 03 C 1 S SA = A6 iC I SA( CM LG ( L ) » CM t G ( L -i ) ,OMEG( L-2 ) , CETT ( L I , DETT (L- 1 ) , 
*UcTT ( L-2 ) ) 

G=C I S 3 A 
D L T = 0 9 

rtRITcI o,250) G, CC , OE T » R CCNT 
102 PRINT i20 

L Me G A=UM LG A*- OE LC M 
IF ( U.MLgA .UT.UMF IN (GO TO 110 
GU Tu .01 
110 wRI TE(6,120) 

12 u FORMAT ( ) 

Q=Q+DfcLU 
OMEGA =UM EGA 1 
I F ( Q. GT® QFIN ) GU TO 1 
Gu TO .00 

C DETERMINE CMLGA FOR EACH u 

300 KUUNT = u 

Uu j.O 1=1,3 

IF(I.Gb®2) QM£GA=OMbGA+CeLCM 

CALL DET ERM{ CMEGA, Q , LM »c I ,NDS, DET , C , EK, W , EQ ) 

xd i=dlt 

Y ( I ) =UME G A 
310 CONTINUE 
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320 C I S SA =28 SC ISAl Y < 3 > » Y 1 2 ) , Y Q ) ,X(3) ,X( 2) »X(I) ) 
0McGa=C 1 SSA 
IHKJuiyT.2w. 0I X0=X(i) 

KQUnT =KUUNT +1 

CALL JET tfRM t OMEGA, W, LM,e I * NpS, GET ,C , EK , W , £Q ) 

X(I) = DcT 
Y ( 1 ) = UME GA 

I F ( KU UN t . GE. 10 ) GO TO 7 u 
i F ( AOS (X 1 i ) ) / A bS ( XO ). LE . A »fc— 5 ) GO TO 34C 
UEL0rt=Dfc LOM/IO . 

CU Ej J 1=2*3 
CMcGA =Q,MEGA + UELCM 

CALL (JET FRM ( GH EGA * U * EM » c I * NOS » DET,C , EK , W * EQ ) 

X { I ) = UtT 
Y l I »= OMEGA 
330 CONTINUE 
GO TJ 320 

34-0 C=0Mi:.tjA*cL**2 
Ww=a*EL*#3*2. 

PRIiyT 250, 0 » U W , GET , KCCNT 
vm R I T_ (b, 23u ) 0, WW»DtT , KCOM 
w = W + JEL 0 

IFIW.G7.WFIN) GO TC i. 

GU TJ a 0 0 

250 FORMA T(/jt2J. 8 ,12x15) 

260 FURMAT ( 1a£9H1u ITERATIONS, NO CONVERGENCE) 

ENO 
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SUbRGUT I NE: DeT ERM{ CMfcGA , Q , EM ,c l ,N OS , OET, C , EK » W,EQ I 

l VaL UAT £ DETERMINANT FCR CumEGA.QI 

DIMENSION EMlil.dlD.wU), E ( Z , 2 ) , C< 1 I , EK< 1 1 

DIMENSION EUIil 

> I u£ Plate EQUATIONS in difference FORM 

W XX I X ) = W ( K + 1 ) -2 . ** ( R I t-fc l N- 1 ) 

X M ( K ) - fc I ( K) #W XX( K ) 

XE M( N ) =XM( K + i ) -2 • * XM l K ) + XM (N-j. i + Q#EQ ( K ) #( W ( K+l I — W ( K— 1 I ) 

-C(X + i )* (W(K + i)-w(N) )+C(lO*(WlK ) - W C K — i I I 
+ ( EK ( K )— CMEGA* 1 « c 2*tM( K) K ) 

ndss=nds +1 

NOit=NPS *2 
NUr.r = NDS+4 
DU ou 1=1, NO E t 

MARCHING PROCEDURE 

M I ) =u. 0 
DO 90 N = 3 » 4 
MN=U 
to( 3 ) =0 • 

W 1 4 I = 0 • 

W<N)=_. 

DU 70 J=3,NUS 

£XM=l. *AMI J) -XMl J-i)-w*cvi(j) *( w( J + il-R ( J~l) ) +C( J + l ) 
i *lw(JH)-W(jn-tlJI*ltilJ)-«tJ-ll )-( FK ( J )-CM£GA**2*EM( J » ) *W ( J J 
tXX= fcxM/t I( J+l) 

Vi( J+2)=£XX+2.*A( j + i i-w(J) 

C ON T I N Ufc 

CUtF FI C I tNTS CF LAST T aU cwUAT IONS 

DU 8u J=NOSS»NDSfc 
MN=MN+1 

t (MN ,N~2 ) =XEM( J ) 

CUNT I NUc 
CUNTINUfc 


GcTERM 1 NAN T OF COEFFICIENTS 

Tl=E(i,l)-.£( I, II 
T2=t:U,E)-E..Ii f 2) 

Ti=T„*c( i j 2 I +T 2*E ( 2 » I ) 

DLT = Ti*T 2 + T3 

Rl T URN 

END 
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Figure 3.- Typical curves from which (flutter) values of q are determined at frequency coalescence. 

(The symbol X denotes flutter.) 
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Figure 5 .- Effect of support stiffness on the flutter mode of the special panel. 
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Figure 6.- Flutter of panels on equally spaced flexible supports. 
(Asymptotic values were taken from ref. 2.) 
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Figure 7.- Flutter of panels on unequally spaced supports. 
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